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ABSTRACT 


This  investigation  is  directed  to  the  analysis  of  a  composite 
disc  in  plane  elastic  equilibrium  The  disc  consists  of  concentric 
regions,  each  of  homogeneous  properties,  end  is  loaded  on  its  outer 
surface  by  an  arbitrary  load.  Reactions  at  the  center  are  prescribed 
to  provide  the  necessary  equilibrating  force  and  moment.  The  equations 
of  equilibrium  are  specific'  in  terns  of  the  radial  and  tangential 
displacements,  giving  a  pair  of  partial  differential  in  plane  polar 
coordinates. 

A  finite  Fourier  transform  is  applied  to  reduce  the  partial 
differential  equations  to  a  pair  of  ordinary  differential  equations 
with  the  radial  coordinate  as  independent  variables.  Inversion  of  the 
transformed  solution  yields  an  infinite  series.  The  analysis  is 
continued  in  further  detail  for  uniformly  distributed  and  concentrated 
normal  and  shear  loads. 

The  deformations  associated  with  the  uniformly  distributed 
normal  load  are  used  as  the  basis  of  a  postulated  friction  phenomenon 
called  the  deformation  coefficient  of  friction.  Assumptions  about  the 
behavior  of  the  strain  energy  as  a  load  is  translated  an  incremental 
distance  give  rise  to  a  friction  force,  vhich  in  turn  permits  a 
friction  coefficient  to  be  calculated.  This  coefficient  is  charac¬ 
terized  by  being  proportional  to  the  load  rather  than  constant.  The 
contribution  of  the  deformation  coefficient  to  the  total  friction 
behavior  is  ordinarily  quite  small. 


vii 
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PART  I. 

INTRODUCTION  AND  HISTORICAL  REVIEW 

A.  Introduction 

This  thesis  is  concerned  with  a  problem  in  the  plane  theory  of 
elasticity.  We  shall  consider  a  composite  disc  consisting  of  two  per- 
fectly  bonded,  coplaner,  and  concentric  regions  which  are  each  isotropic 
and  homogeneous.  The  di3c  is  loaded  by  arbitrary  stresses  on  tne  outer 
boundary.  A  concentrated  central  force  and  moment  equilibrate  the 
boundary  tractions.  A  formal  solution  in  the  form  of  an  infinite  series 
i3  obtained  for  the  general  case.  Solutions  are  obtained  explicitly  for 
a  few  special  loadings,  and  numerical  values  obtained  for  the  particular 
case  of  a  normal  stress  of  unit  intensity  distributed  over  a  segment  of 
the  boundary. 

A  friction  mechanism  is  postulated,  based  on  the  strain  energy 
stored  in  the  system  during  an  incremental  translation  of  the  area  of 
application  of  the  normal  stress.  This  deformation  coefficient  of 
friction  depends  on  the  elastic  behavior  of  a  smooth  disc  under  the 
influence  of  a  normal  load.  It  does  not  take  into  consideration  other 
possible  sources  of  friction,  such  as  sticking  and  the  interaction  of 
asperities.  The  deformation  coefficient  is  ordinarily  quite  small 
compared  to  Columb  friction,  and  has  the  unusual  property  of  being 
proportional  to  the  load.  The  frictional  force  contributed  by  this 
mechanism,  therefore,  varies  as  the  square  of  the  normal  load. 

The  mathematical  analysis  used  is  somewhat  unusual.  The 
equilibrium  equations  are  established  in  terms  of  displacements.  A 


>  • 
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finite  Fourier  exponential  transform  is  used  to  reduce  the  system  of 
partial  differential  equations  to  ordinary  differential  equations.  The 
transformed  solution  is  inverted  by  an  infinite  series.  The  transform 
method  has  properties  which  make  it  quite  useful  in  circular  regions. 
Unlike  the  more  familiar  Kolos  of  f-Muskhelishvili  method,  in  which 
complex  potentials  of  known  form  are  required,  the  transform  method 
leads  to  the  solution  through  a  direct  process.  The  resulting  infinite 
series  corresponds  to  the  series  required  to  establish  the  harmonic 
functions  in  a  composite  region  when  the  Kolosoff-Muskhelishvili  method 
is  used.  In  the  particular  case  when  the  disc  is  homogeneous  the  series 
solutions  may  be  exhibited  in  closed  form.  When  the  distributed  load  is 
condensed  to  a  concentrated  load  the  solution  agrees  with  previously 
known  solutions. 


B.  Historical  Review 

The  analysis  of  circular  regions  in  a  state  of  plane  elastic 

equilibrium  has  been  the  subject  of  substantial  investigation.  One  of 

22* 

the  earliest  reports  on  the  subject  was  that  of  Michell  in  which  he 
examined  a  class  of  problems  in  the  plane  disc.  Michell  considered 
concentrated  load3  in  the  interior  of  an  infinite  plane  and  on  the  edge 
of  a  plane  half-space.  By  superposition  of  these  effects  he  developed 
formulas  for  the  disc  loaded  at  the  center  and  edge  by  concentrated 
loads.  He  gave  explicit  formulas  for  the  stress  field  due  to  concen¬ 
trated  loads  at  the  edge  with  an  equilibrating  reaction  at  the  center. 


•Throughout  this  thesis,  superscript  numbers  refer  to  the 
similarly  numbered  items  in  PART  VII,  BIBLIOGRAPHY. 
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These  formulas  are  unfortunately  incorrect  due  to  failure  to  properly 
annihilate  the  reactions  on  the  boundary  due  to  the  edge  load. 

Mindlin^  extended  Michell's  problem  to  consider  the  effect 
of  translating  the  central  force  to  an  arbitrary  interior  point, 
directed  along  a  diameter.  As  a  starting  point  he  used  Michell's 
problem,  for  which  he  expresses,  without  derivation,  the  correct  Airy 
stress  function. 

The  development  of  the  Kolo3off-Muskhelishvili2^*  complex 

variable  method  made  it  possible  to  solve  a  large  class  of  problems  in 

the  disc.  In  most  cases  it  is  necessary  to  use  a  catalog  of  harmonic 

functions  to  satisfy  the  boundary  conditions.  Mikhlin*"’  examined  a 

composite  disc  loaded  on  the  boundary  by  self-equilibrating  forces.  His 

solution  is  in  the  form  of  infinite  series  of  harmonic  functions, 

q 

Dokos  examined  a  composite  consisting  of  a  disc  loaded  at  the 
center  and  embedded  in  an  infinite  plate.,  His  method  consisted  of  super¬ 
position  of  field  equations  for  stress  and  displacement.  Dunaurs  and 
Hetenyi11'  ^  '  10  examined  a  series  of  problems  of  the  elastic  disc 
embedded  in  an  infinite  plate.  The  Kolos of f-Muskhelishvili  formulas 
were  used  somewhat  in  their  development,  but  the  method  was  basically  a 
superposition  of  appropriate  Airy  stress  functions. 

Ng^®  resolved  Michell's  problem  using  the  complex  variable 
method.  His  results,  given  in  terms  of  displacements,  can  be  shown  to 
be  in  agreement  with  those  of  Mindlin. 

20 

A  recent  investigation  by  Maye  and  Ling  used  finite  Fourier 
transforms,  along  with  the  theory  of  couple  stresses,  to  examine  the  ring. 
Maye  developed  the  equations  for  a  composite  disc  as  a  special  case. 


k 


Due  to  the  complexities  associated  with  couple  stresses  an  explicit 
analytical  solution  vas  not  given.  His  technique  of  inversion  uses 
only  the  real  part  of  the  series. 

8 

Finite  Fourier  transforms  have  been  developed  by  Doetsch 

1  A 

and  KniessA  •  Applications  of  these  transforms  to  various  vibration 

2  3  1* 

and  heat  conduction  problems  have  been  developed  by  Brown  *  *  and 

Roettinger^*  Generally  these  investigations  have  been  restricted 
to  the  sine  or  cosine  transform. 
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PART  II. 
THEORY 


A.  Formulation  of  the  Problem 

Consider  a  circular  region  in  a  state  of  plane  elastic  equilib¬ 
rium.  The  equations  of  equilibrium  in  cylindrical  coordinates  may  be 
written  in  terms  of  the  radial  and  circumferential  displacements 


(A*  V 
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For  a  composite,  as  in  this  analysis,  Eqs.  (l)  must  hold  in  both  regions. 
Constants  and  variables  associated  with  the  core  will  be  distinguished  by 
a  bar  (e.g.  U  ),  while  the  unbarred  symbols  will  apply  to  the  ring,  or  to 
the  whole  system,  as  appropriate. 

The  elastic  constant  It  is  related  to  the  more  familiar  Lame'’ 
constants  and  Poisson's  ratio,  in  the  case  of  plane  strain,  by 


Jk=(A+f)/p  =l/(l-2v) 


(2) 
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and  for  plane  stress  by 

J,  =(3X+2r)/^  =  (l+M)/a--*)  (3) 

Stresses  in  each  region  may  be  expressed  in  terms  of  the 
displacements  by 


cr*  = 


(M 


We  shall  take  as  our  dimensions  a  disc  whose  outer  radius  is 
unity  and  for  which  the  radius  of  the  interface  is  8.  Prescribed 
surface  tractions  on  the  outer  circle  are  equilibrated  by  concentrated 
force  and  moment  at  the  center.  The  polar  coordinates  (f,-©)  of  a 
generic  point  are  measured  from  the  center  of  the  disc  and  the  X  axis. 
Figure  1  is  a  representation  of  the  elastic  system. 


Figure  1.  Composite  Elastic  Disc 
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It  is  convenient  to  separate  the  boundary  stress  into  its 
normal  and  tangential  components.  Therefore  we  express  the  conditions 
at  r  =  1  by 


ov  ( 1,0)  =g  0Q-) 

r*.  Ci  ,$)  =  -f  600 


(5) 


In  order  to  establish  tne  conditions  at  the  center  of  tne  disc 
we  consider  the  retirement  that  the  concentrated  central  force  and 
moment  equilibrate  the  force  and  moment  resultant  of  the  boundary 
stresses.  The  force  will  be  compactly  expressed  in  complex  form  terms 
of  its  X  and  y  components.  The  force  and  moment  resultants  are  then 
given  by 


fTT 

Fx  *  *  =  J  6  dL-0- 

li  =  \n  ffe)  dO- 

IT 


rfe  will  first  look  at  a  central  hole  of  radius  £.  Consider  the  surface 
of  the  hole  to  be  acted  on  by  a  constant  shearing  stress  and  by  a 

normal  stress  The  boundary  conditions  associated  with  these 

requirements  are  given  by 


(6) 


o-r(€.,-90  =  g  (£) /e 

(7) 

r^(e,h)  = 


B 


The  condition  of  equilibrium  with  the  surface  tractions  requires  that 

Z  §C<&  =  r.  + 1  f, 

(8) 

y,  f  <&  =m  =>  F  =  (1/2TT 

The  condition  of  concentrated  force  and  moment  at  the  center  is  obtained 
from  the  limit  as  O  . 

The  boundary  problem  is  completely  prescribed  when  we  add  the 
conditions  that  U,  nT,  ,  and  are  continuous  at  the  interface 

between  the  regions. 

Eqs.  (l),  (M ,  (5)»  and  the  limit  of  (7)»  subject  to  the  side 
condition  (8)  and  the  matching  conditions  at  the  interface,  establish  a 
boundary  value  problem  for  the  disc.  We  shall  examine  the  problem  in 
the  general  case  and  then  look  at  some  particular  examples  in  detail. 

fl.  Transformation  of  the  Problem 
The  solution  of  the  posed  boundary  value  problem  will  be 
obtained  by  means  of  the  finite  Fourier  transform.  Our  definition, 

g 

which  is  similar  to  that  of  Doetsch  ,  is 

2.  =  =TW  C  Z<W  e"**  dG-  «> 


whose  inverse  is 


9 


OP 

zceo  =£zs  e 


(.S'© 


(10) 


S=-«» 


The  correspondence  of  Eqs.  (9)  and  (10)  is  well  known  from  the  theory  of 
7 

Fourier  series  .  Except  for  the  fundamental  analysis  by  Doetsch  and  the 

20 

recent  work  of  Maye  ,  however,  the  use  of  Eq,  (9)  as  a  transform  is 

rare.  The  corresponding  Fourier  sine  and  cosine  transforms  have  been 

used  somewhat  more  frequently  for  the  solution  of  initial  and  boundary 
2  3  1*  30 

value  problems  ’  *  *  .To  obtain  the  transforms  of  derivatives  we 
apply  the  definition  to  the  T\  th  derivative  of  z&).  Integrate  once 
by  parts,  to  obtain 


[zfr°  ©■))  =  Ls  (z^  0(9-)] 

[Z'-’CTT)  -  Z(*-°(-TT)] 


(11) 


in  a  circular  region  ■©"  =  IT  and  ©■  =  ”7f  coincide.  Tnerefore  if  ZC-9-) 
and  its  first  TV-1  derivatives  are  continuous,  the  constant  part  of 
Eq.  (ll)  vanishes.  Transforms  of  any  derivative  of  z(e-)  can  then  be 
written  in  terms  of  Z s  as 


=  (is)*  {Zflw}  =(Ls)hZ: 


(12) 


The  power  of  the  finite  Fourier  transforms  in  circular  regions  is 
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enhanced  by  this  simplification.  For  any  function  which  is  reasonably 
veil  behaved  it  is  possible  to  orient  the  coordinates  to  assure  the 
required  continuity.  In  the  analysis  which  follows  it  is  assumed  that 
Eq.  (12)  holds. 

The  transformed  equilibrium  conditions  Eqs.  (l)  are 


C*  +0  <  +  Uj 

+  i  =0 


-  i.  //+  -i-  rr'  -  ^4'^f  ~l~  *  nrj  =  0 

This  coupled  pair  of  ordinary  differential  equations  may  be 
solved  by  elementary  methods,  giving  for  |S  | &  2 

»♦!  .  i  -(*+!>  .  .  ,  «-»  .  ,  ,  -(4-1) 


•Us  =  U„  r 


+  UM  +  UJS  r*’1  +  Us4  r 


;  [^5+2*  ■♦■2,1  r  SM  II  P 

^  =,L  u«  r  ”u“r 


-(S  +  1) 


(lM 


~  U44 


+  U„r 


4l»+2 
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For  S  -  0  the  equations  become  uncoupled, 

uo/  +T  uo/-'F  uo  =0 

(15) 

<'  +T-^--F«ir0  =0 

vnose  solution  is 

^0  “  ^01  f  +  ^oi  r 

(16) 

^0  ~  ^01  r  4  Xu  r 

Wnen  |S|=1  the  solution  £iven  by  £qs.  (lU)  is  incomplete 
because  of  a  repeated  function.  However,  if  we  rewrite  bqs.  ( 13 )  for 

S  =  ±  1 

O  +  Ou/^  +is[i-nrs/- i^nrj  =0 

(IT) 

V4-  r  <  -  <vj  =0 

we  ootain  the  solution 


*  n  /  i  U2. 
r  *  r1 


U.  -  is 


\A  s  -  U41  r  Z  +  Usi  r  1  +  U53  +  Us^ 


at5-  is 


U  r*  -  U  r’*  +  (J 

•A-Z  Si  52  S3 


+  R^+Us^r] 


(18) 
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The  transformed  stresses  are  expressed  in  terms  of  the 
transformed  displacements 


f  1 

[(>  +  1)  u/  +  r  (U4+LShr4)] 

* 

(0 

If 

(A-l)u/  +  -^-(us  +  is  ^S] 

Ts  -  ~~  T“(SUs  +  i.AT$)J 

More  explicitly,  for  5  =  0 

crro  =  2y*aU#1  -  UM  r'*) 
a„  =2/aUol  +  Uolr-2)  (20) 

T0  —  -2^»  Y01  r 

tfhen  S  =±  1 

^  =  -*>  [x^Um  r  +  2  U«  r'S  r'1] 

=-2r[jrrUsl  r  -2UMf 3  +PJ-LU  r  ] 

TS  =2isr[l^^r 
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For  $  £  2 
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>S-2  US1 


s 
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-(S-M)U 


si 


-(5  +  2) 
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+  (s-i)U53 
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o;5  =  -2/' 
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(22) 


+  (S-1)U  r*'1  -  y  r-S 
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ir 
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C.  Solution  and  Inversion 

The  solution  of  the  disc  problem  requires  tne  determination 
of  eight  constants  for  each  value  of  5  ,  four  in  each  of  the  ring  and 
the  central  core.  The  eight  boundary  conditions  to  be  satisfied  are 


^  CO  -  9S  =  Tri  L  9»)e"‘4d^ 
TgCO  =  =^f  1„  f(£) 


(23) 


1U 


<Trj  (i)  -  0"rj  (%)  =  0  (c) 

Ts  M  =  0  (d) 


■Us(0  -  us(0  =  0  (e) 

(23) 

ATs(S)  ~  ATS  Cf>)  =  0  (f)C°nt 


where  (g)  and  (h)  are  to  be  satisfied  in  the  limit.  Evaluating  the 
integral  in  (h)  gives 


(s  =  o) 
(s  =0) 


(2b) 


We  will  examine  the  boundary  conditions  for  three  classes  of 
values  of  S,  that  is,  for  S  =  Ot  |S|  =  1  ,  and  |s|^2.  • 

If  we  express  the  boundary  conditions  explicitly  for  S-0 
we  can  write  for  the  Y0^  constants 


V0x  =  -f  oA/*  <b) 

(25) 

2  ^02.  =  ^ 


(d) 
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+ 

s"  vw  -  SV01  -  S"  V#l  =  0 

(f) 

(25) 

Cont 

e'1  v 

c  voz 

=  Ke'V2? 

(h) 

where 

=  f0  =  n/2TT 

(26) 

Conditions  (b)  and  (h)  arc  redundant.  As  a  result,  one  constant  (we 
choose  V01  )  oust  remain  undetermined.  The  undetermined  constant  corres¬ 
ponds  to  a  riGid  rotation  of  the  system.  The  constants  are  found  to  be 


VQ1  =  -  M/4TT)u 

Vol  —  -M  A-TT^:  (27) 

V01=nCyu tY01 

The  boundary  conditions  for  the  UQj  coefficients  are  given 
explicitly  by 

^^01  “  ^02.  ~  9o/2^  (a) 

(28) 

^01  ^02  “  ^01  +  F  k  ^01  =(^  (c) 
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su01  +s"uM  -  f>u0l  -$>'‘U0l  =  0 

^01  ^  ^02  ”  9o  ^ 

Application  of  the  limiting  condition  to  (g)  shows  that 


The  remaining  conditions,  expressed  in  matrix  form  are 


1 

o 

?- 

1 _ 

^01 

'  r 

ywl  *  ~/U  & 

cn)0^ 

II 

0 

.S*  1  -s1 

A. 

.0. 

from  which  we  find 


_  9oO  *pX') 
uoi  a /.Do 

_  % 

uoi  ~  zrD0 


iD0 


(e) 

(28) 

Cont. 

(6) 


(29) 


(30) 


(31) 


where 


d0  =/.*(  i  -s1)  i  au1) 
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Consider  next  tne  conditions  :'or  Isl  =  1  .  dubsti tut  in j 
nqs,  (16)  and  {21)  into  oqs.  (23 )  *ives  tne  explicit  conditions 


2‘  u..  -  ?-  uM  -  25  +  1 


J k-2  SI 
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(  c) 


1*1 


$4 


=  0 


-1 


,  +2/.5'*U,,  +  d^-Ui4 

-ft-2  s*  r  il  Jk*i  44 


-3  rr 


_  iMAu  .  -2pi"U.,  -- £ 

1-2  si  ”  44  A 


A  S4 


=  0 


(32) 


(d) 


S  Usi  +  S  Us2.  +  US3  +  S')  Us^ 

-  ^  u51  -  s  US2  -  U53  0  us5  -  o 


(c  ) 


s'1  U„  +u 


si 


SI  x 


°S4  + 


5*4 


(?) 
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-2±L  u 

i  -2 


*+2.  US4 


2JU 

TT 


From  Eqs. 


+  2«'iUsi 

(6),  (8),  and  (9)  we  find 


0 


(e) 

(32) 

L- ont  • 

(h) 


3s  =  9s  -  isfs 

Combining  Eqs.  (32a)  and  (b)  gives  \J&+  ;  (g)  and  (h)  give  Uj,,  . 
Taking  (g)  or  (h)  in  the  limit  as  £ —  0  gives  Ui2  .  We  obtain, 
thereby 


U„  =0 

,,  .  g+iK3,-lSf.) 

U34  “  4^a+l)  (3M 

y  _  (i^XV^s) 

i4  4^0* +1) 

Elimination  of  these  constants  from  Eqs.  (32)  leaves  five 
equations  in  five  unknown  constants.  The  constants  (J53  and  (J53 
appear  only  in  combination.  We  can  eliminate  their  difference  by 
combining  (e)  and  (f).  The  remaining  four  equations  have  three  unknown 
quantities.  The  equations,  however,  are  redundant.  One  may  therefore 
be  eliminated.  We  may  write  the  remaining  three  equations  in  comoact 


form.  Let  us  define 
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X51  - 


4U 


SI 


Jk-2 


X52  = 


0 


51 


9s 

i^Ui-iXi+D 


S 


+1) 


(35) 


zr,  =1  ^  =  -(2A+l) 


.Ve  mac,'  now  express  tne  remaining  equations  in  matrix  form. 


A  1 

0 

x.; 

2 

'  ar.^cJE*0 

y  M  S  yU 

•pi  i"* 

*sx 

■2»« 

8  yyi  (A  -  -4  ) 

Cfc  +  2)S4  -1 

-d  *1)1* 

*53. 

i-1 

S  [)*  Jk(i+l)-y*A(fc«-lJ 

The  solution  to  Eqs.  (36)  is 

US1  = 

1  r-i 

+  +Jd*  )]Qsj 

os,  =  It/-*  t2i_ 

1  i«‘ 

}Qsj 
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_  ya+ntt-a) 

4Dtgl 


yild-s4) 


-2yi(J  +  l)il-ai+lI -2i)s4]  jos. 


(37) 
Cont . 


where 


D,  =>.Jtd+2)(l-S4)  +fl6fe+JkS‘*+2i4) 

One  constant  remains  arbitrary.  (We  choose  u  53  )  •  The  remaining 
constant  is  then  given  by 

U4J  =  U  S3  +fusi-u41')s1  +  U51fc'*  +(U^-Ust)JUji  (36) 

We  have  expressed  the  solution  for  S  =  1  1  .  The  two  terms 
in  the  solution  for  each  of  the  constants  correspond  to  the  normal  and 
shear  components  of  the  boundary  tractions. 

Let  us  next  examine  the  constants  for  |<S|  £  2. ,  We  observe 
first  that  if  we  use  £qs.  (22)  to  explicitly  state  boundary  conditions 
given  by  Eqs.  (23g)  and  (h)(  we  obtain  relations  of  the  form 


n  -  ^  !<*♦*>  n  i*  T-7  _  s*i  — 

US1  -  Qte  U»i  +Qxe  Us,  +a,£  q 


S-l  — 


0*  =  b,clS  u.t  +  b,e,,”uUSI  +  bj£‘"‘g  i 

where  and  bj  represent  collections  of  constants, 
limit  as  £— VO  for  5^  2.  ,  we  find 


(39) 


uSJ  =  us4  =0 


(ko) 


If  we  take  the 
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Corresponding  statements  can  be  written  for  $$  -2  ,  from 
which  we  find 

DSJ  =  US1  =0  (si-2)  in) 

As  in  the  case  of  the  indices  S  =  0  and  5=1  1  ,  it  does  not  matter 
how  the  normal  stress  is  distributed  at  the  center. 

For  ease  in  handling  the  equations,  it  will  be  convenient  to 
express  the  indices  in  terms  of  the  positive  integer  t.  *  defined 

t  =  |s|  <«) 

We  can  now  take  advantage  of  the  symmetry  of  form  for  positive  and 
negative  indices  by  defining 


[xtl 

+l) 

[U*. 

U.44] 

[X„ 

x.»] 

=  ^ft+1) 

[U  t» 

^-ta] 

[Xts 

X.„] 

=  r(t-l) 

['l.J 

^-ii] 

[*.4 

X-44] 

<£  + 
cJK  +* 

II 

[*•*** 

U-ti] 

X-tj] 

=  Tw 

[D«. 

D-tJ 

[*« 

X-tiJ 

=  f(t-i) 

[0*3 

u.„] 

(»*3) 


We  also  define 
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=  S 

QSi  =  -  L  (s/t)  fyA 
A  j  —  ~  1  A  2  —  1 

If  we  recombine  the  remaining  six  conditions  of  Eqs,  (23)  in  pairs  ve 
may  write  the  boundary  conditions  explicitly  and  compactly  in  matrix 
form,  as  is  shown  in  riq.  (1+5). 


( 1+1+) 
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2k 


The  solution  ol’  hq.  (U5)  may  be  written  in  relatively  compact 
form  if  we  define  the  following  collections  of  terms: 

w  =pIU+2)  ->JKI+2) 

X  —  yl  -  j*  A 
Y  -Y.? 

Z  = 

Z  =  JO.  (1+2)  +j:l  <1" 

Dt  =  -HZ  +  tlJk‘YZ  hUt'° 

-2(t‘-t)l2YZ  &** 

+[tiiirz+*/(Mi)x]r\AYwr 

The  solution,  then,  is 
Z 

Xsl  =  45^  (M2*  -(t+l+AjJA'YZ 

j;1  ,  —  3  t  ) 

+(ttl)A  YZ  8  ] 

X»  ='f2Qsi  (t-O- 

-  Aj  [t’A!  YZ  +  ‘kr(JUi)x]sl 

-ct+H-AjUrwV*) 
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x«  = -572°s^[Aja'°  'I]AZZ 

i=l 


+  [i‘*‘YZ  +  4 f  (Jt  +  f) *]  s"4*”] 

xS4  -  -^2°^  f  [AJ(t  ■ 0  -*]  Yf 

j=t 

-Ajd-OJt’YZ  S1  -iYWSU“l>l 


X 


55 


ifiia up 

»t 


,  2r*-o  ,*■) 

-ft  ♦UAjJiYJ  +(tM)JkXS  ] 


(W 

Cont. 


-2Aj('t-i)Xi1  -ft-lXt  +1  +Aj)JUl  Yi" 


r  +  2x]  s1^1’] 
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The  solution  is  now  formally  complete.  The  coefficients  have 
all  been  calculated,  and  may  be  summed  to  give  the  desired  functions. 

The  inversion  formula,  Lq.  (10),  is  somewhat  more  manageable  if  broken 
up  into  sine  and  cosine  series.  Inversion  of  the  results  will  therefore 
be  effected  using 

(1*8) 

ZO.e)  =  Z0(V) 

oo 

+}_  { [Z,(D  +  Z.s(r]  C0*ao)  +  t[Zs(r)-Z.s(r)]  wfs«)j 

S-  1 


HAST  III 


uxa:lplej 


A.  Uniform  and  Concentrated  formal  Loads 
As  tne  first  example  we  shall  look  at  the  composite  disc 
loaded  in  its  plane  by  a  unit  normal  stress  over  the  arc  (-<*>/ 2, 4>/2  )  . 
The  boundary  conditions  may  be  explicitly  expressed 

(l  (-■£-<£<-?-) 

crf(i,e)=  g  (el  =  j 

10  elsewhere  ( i»9 ) 

T^Cl.G)  =  {($)  =  0 

The  corresponding  conditions  for  the  transformed  problem  are 


f  <p/zjt  (5=0) 

(54>/a)/TTs  (s  *  0) 

f*  -- 0 


i\  li  mi  tine;  case  of  tne  distributed  normal  load  is  tne  con¬ 
centrated  unit  normal  force.  To  obtain  this  case  from  the  distributed 
load,  we  first  diviac  by  2 .44/rv  to  normalize  the  force. 

As  4»  vanishes  the  load  approaches  a  concentrated  force. 
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The  boundary  conditions  for  the  natural  problem  can  be  written 


($-)  =  { 


cp-0 


l  o 


elsewhere 


(51) 


The  transformed  condition  for  the  concentrated  normal  force  is  therefore 


S  s  ~  2  TT 


fs  =  o 


(52) 


Substitution  of  these  values  into  Eqs.  (29)  and  (31)  gives, 
for  5-0 

Distributed  Load  Concentrated  Load 


(Joi  - 


4-V  D0 


u.i  = 


+ 1) 
4-TTD0 


u.,= 


I)  -  II  = 

U02  '  4-TT,.  U01 


+  TT^D0 
4TTrD0 


0  =  -A±i 

u0i 


4TT  Dc 


(53) 


u02  =  0 


U0J  =0 


rt'e  will  assign  the  value  Vjjj  =  o  to  the  arbitrary  constant.  Since 
the  moment  M  is  zero,  Eqs,  (27)  give 


v.  -  V.  =  o 

Yoj  Toj 


(5*0 


The  constants  required  to  evaluate  the  coefficients  for 
| $|  =  1  .xe  obtained  from  Eqs.  (3M,  (37)*  and  (38).  Eqs.  (50)  and 
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(52)  show  that  for  botn  the  distributed  and  the  concentrated  normal 
load 


9-5  "9s  f*  =  0 

The  constants  needed  for  5  si  1  are  obtained  using  the  quantities 
defined  in  Eq.  (35)«  in  particular,  we  find  for  the  present  case 


Si 


0-12  -  0,2  =  0 


We  will  define  the  arbitrary  constants 


(56) 


U.I3  =0U  =0  (57) 

Witn  these  results,  we  observe  from  Eqs.  (31*)  and  (37)  that 

u-«i  =  uij  O-tj  =  UtJ 

and  that  these  constants  are  all  real.  From  Eqs.  (18)  and  (21), 
tnerefore,  we  find 

U.,  =  U,  ^r,-i  —  = 

(59) 

OT.j  =  -  nrt  T_i  =  -  T, 

for  the  barred  as  well  as  the  unbarred  values.  We  note  also  that  the 
even  coefficients  are  real  and  the  odd  coefficients  are  imaginary. 


In  a  similar  fashion,  using  Eqs.  (Uo),  (i*3),  (M)  to  find 
the  constants  for  l^j  ^  2  ,  we  find  that 
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U-SJ  =  Us<5.  j)  U.5J-  =  UsCs.p  (t>0) 

In  other  words,  the  constants,  associated  with  corresponding  positive 
and  negative  indices  and  whicn  multiply  the  same  power  of  T  ,  are 
equal.  Therefore,  Eqs.  (ll*)  and  (22)  show  that 

^  -s  ~  ^5  crvs  -  o.r5  - 


Here  also,  the  even  and  odd  coefficients  are  real  and  imaginary, 
respectively. 

The  inversion  of  the  transformed  solution  requires  only  the 
cosine  terms  of  Eq.  (1*8)  for  the  even  coefficients  and  the  sine  terms 
for  the  odd.  Thus 


r\r 


SM 


rvrs  (.S©-) 


(62) 


with  corresponding  formulas  for  the  stresses. 

The  solution  is  now  completely  determined  in  series  form. 
The  necessary  coefficients  are  given  by  Eqs,  (lk),  (16),  (l8),  (19)» 
(20)  and  (21).  The  boundary  values  are  given  by  Eq.  (50)  or  (52) 
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depending  on  tne  load. 


o.  Unii’orai  and  Concentrated  dhear  Loads 
Tne  loading  of  a  disc  by  a  unit  snear  stress  distributed 
over  tne  arc  (  -<*>/*  %<x>U  )  jives  the  boundary  condition 

.  ,  1  (-$  <■&  <$) 

Trft(  1  ,*)  =  f  (e)  = 


0 


elsewhere 


cr/i.O-)  -  g(eO  =  0 


file  boundary  condition  for  the  corresponding  concentrated  load  is 


T  r$(l  ,-9)  =  -f  (-0)  = 


o 


elsewhere 


(63) 


(61*) 


The  transformed  boundary  conditions  are 

2  TT 


tts 


(5  =  0) 
(5*0) 


(65) 


and 


r  i 

+  s  "  TTT 


(66) 


for  tne  distributed  and  concentrated  loads,  respectively.  In  each 


case 


9»  =0 


(67) 


^e  now  can  determine  the  coefficients  needed  to  evaluate  the 
;eries  for  displacements  and  stresses.  If  we  assign  to  the  arbitrary 


t  /. 


01 


const an 


the  value 
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^01  -  ■  ^01 

we  can  write  from  Eqs.  (27) 

V0l  =  V2/* 

V02  =  --fo/2/* 

Vo.  = 

Ki  -  ~  fo/ 2.  /* 


(68) 


(69) 


From  Eqs.  (29)  and  (31)  we  find  that  the  coefficients  U  oj  all  vanish. 

To  find  the  coefficients  for  S  =  Jl  1  ,  we  observe  first  that 
the  only  non-vanishing  right  hand  members  of  Eqs.  (32)  are  in  (b)  and 
(g).  dince  these  are  odd  in  5  and  imaginary  the  constants  U4j 
and  Ujj  will  have  the  same  property,  if  we  again  define 

U,J  =  U_u  =  0  (70) 

From  Eqs.  (l8)  and  (21),  therefore,  we  find  that  ,  C7”r5  ,  and 

are  odd  in  5  and  imaginary  while  and  are  even 

and  real. 

For  |<s|^  2.  ,  a  symmetric  or  concentrated  shear  load  is 
seen  from  Eqs.  (40),  (43),  (44)  and  (47)  to  have  the  property 


=~U3(5-j) 


U-Sj  =  “  ^S(5-j) 


(71) 
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*:ierefore  U.j  ,  0~rS  ,  and  CJ^s  are  °Hd  in  S  and  jure  imaginary 
waiie  hf^  and  Ts  are  even  and  real.  Combining  the  terms  in  the  form 
snow n  by  Hq.  (U8)  we  obtain 


u 


or 


s-i 


nr0  +  Z 


-dfn.  (s&) 


cxro. 


(s4) 


s=i 

wnc-re  tne  coefficients  have  all  been  determined.  Similar  expressions 
no  Id  for  the  stresses  and  the  formulas  hold  for  both  tne  unoarred  and 
oarred  variables. 


(72) 


C.  Remarks  on  an  Arbitrary  Load 
The  examples  just  cited  were  chosen  for  illustrative  purposes 
because  of  their  simplicity.  However,  any  load  may  be  treated  by  the 
same  technique.  He  can  break  an  arbitrary  load  into  four  components, 
that  is,  the  even  and  odd  parts  of  the  normal  and  shear  components.  The 
even  part  of  the  normal  component  and  the  odd  part  of  the  shear  component 
load  to  coefficients  IX  s  ,  CTri  •  and  which  arc  even  in  S 

and  real,  and  to  nTi  and  which  are  odd  and  imaginary.  The  odd 
part  of  the  normal  component  and  the  even  part  of  the  shear  component 
lead  to  coefficients  witn  the  opposite  character.  Thus  each  function 
when  expanded  in  the  form  of  £q,  (U8)  is  explicitly  real.  The  solution 
for  any  load  can  then  be  written  in  series  form,  where  the  coefficients 
are  given  by  the  formulas  developed  above, 

D.  The  Homogeneous  Disc 

Tne  homogeneous  disc  may  be  treated  as  a  limiting  case  of  the 
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conposite  disc  in  several  ways,  sucn  as  letting  tne  interface  radius 
8  vanish  or  become  unity,  or  by  assigning  the  same  set  of  values  to 
the  material  constants  in  each  region.  Alternatively ,  the  Homogeneous 
disc  may  be  treated  as  a  boundary  value  problem  in  one  region.  The 
technique  is  straightforward  by  any  of  the  methods.  The  solution  has 
been  obtained  both  by  solving  the  problem  directly  and  by  equating  the 
material  constants.  The  agreement  of  the  so.wtions  provides  a  check  on 
tne  calculations,  because  of  the  relative  simplicity  of  the  series  co¬ 
efficients,  we  will  write  then  explicitly,  de  prescribe  the  values 

U-u  =^.>  =  0  V01  -  -  Y01  (73) 


to  the  rigid  displacements,  tie  obtain  the  values  for  tne  coefficients  as 
follows  s 
For  S  =  0 

u0  =gor/V  AT0  =  -f0(l-r’) 

®Yo  0o  T0  =  f0/r 

0*0  ■  9° 

For  5  -  ±  1 

^s-5ron5[«-2)rl-z-*fl+2)1^r] 

+2*a*2)-M 


(75) 
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.  9s  r  f. , xl>,r-<i  ,  ojumsf,  f  -iv 

-  2a+o Lr i'^2i+1^r  J  2a+n  (r "  r  > 


^o-r-*)  *  rrrtn  b(u +ih  +  f‘] 


(75) 
(Cont. ) 


Ts  --173TTrfr-r")  *7c^ir[^  +  l)r  +r-] 


For  5  2  2  it  will 


be  convenient  to  work  only  with  positive  values 


of  5  .  */c  combine  the  positive  and  negative  values  of  the  index  by 


defining 


Gs  =  9*  +  **9-1 


fs  =  f,  +  rf  ( 


G5=gs-3-g.s  = 


where  the  symbol  3*  is  assigned  the  value  plus  or  minus  one,  as  appro¬ 
priate.  The  coefficients  for  the  series  are  then 

u>  +aru->  =  - Tpfem  +  '^s)rsM 


[&kt  +  i(s-2)F$]r5'1 


aj\  +  3*  aT_s  = 


•.TfcjGG,  -te-2)Fjr- 
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crrs  ♦  fcr,.,  -  -£-[-(s-2)(Gs  +  i  F,  )r‘ 

+  [sGs  +  i(s-2)  Fs]  r* 
<T*S  +  <r cr^.t  -  -J-((s+2')  (G,  +  i  Fs  )r‘ 

- [s&,  +  Us-2) F,]rs  *  ] 
t4  ♦  iTT.,  =  ^  {-s  (iGs  -  Fj)  r s 

+  [is  G,  -(S-2)  Fj  r*’*  ] 


We  have  now  exhibited  explicitly  the  coefficients  required  for 
the  solution  in  series  form  for  any  loading  on  a  homogeneous  disc. 

For  the  four  specific  loads  described  above,  we  can  carry  the 
solution  a  bit  further.  For  the  unit  normal  load,  distributed  over  am 
arc  ,  the  explicit  solution  is 


_i _ r  **■ 

U  -  4  IT  \  T 


-[a-2k4  +  2xa+2)j^r] 

r  i_1]  Ocm-Gs^)] 

sTz 


(78) 
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+  2£E4^fir"‘  -  W  r‘']^(^/2^^)j 

c  -  •% 


CTr=i{f  ^[rt(2A,0r-^y^ 


&*-2 


r*  )uum,(s4/2)  co^C^)  ] 


r  —  1  ^  ^ a  f  y.  ^  4vr*  (4>/i)c<re^ 

rrl  2  '■,1  _  '  *  +  l 


♦  7  (^  r  *  —  rs  *  )  ~&4sn,  (b^/l)  Co n(s&) 


(78) 

( Cont . ) 


Tr#  =1V((r-r~1)‘%y)^4 


S:2 


(r  —  r  )  *&Lnx  -4-t/rv  (sQ") 
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Let  us  define 

oc  =  4^/2  +  O  (5  =  <P/, 2  —  0  (79 ) 


Combining  the  terms  of  £\,  (77),  we  obtain  series  expressions  of  the  form 


C(i\  go) 
S(r,  go) 
K  (r,  co) 


rs  Cc^(sco) 
s 


>*  -^Ltv  (so) 

S 


f  ^  C/JQ>  (sGo) 


(80) 


2(r,co)  -  y  r 5 ( sco) 

s=i 


These  series  are  all  convergent  inside  the  unit  circle,  to  the  formulas 

17 

given  by  Jolley  . 

C(r,Co)  —  “ (l/2)^laj  NCTjCo) 


5(  r,  go) 

K(r,co) 

ZCr.o) 


jbdsn. 


-i 


t4Wv  CO 
1  -  r  Ora.  as 


t  ( CXrQ'  CO  -  O 
H(r,oo) 


r  ^tv  go 

H(r,co) 


(81) 
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wnere 


M  (r  ,co )  -  1  4-  r2  -2rccraco 

In  the  sane  fashion  we  can  examine  each  of  the  other  loadings, 
lor  the  concentrated  loads  we  need  two  more  series  expressions. 

Q(r.CJ)  =^>  srscWso>)  =  r 

S'-\ 

(82) 

P(r.co-)  =>  sr$^wi(sco)  =  ■rQ-ra)Wco 

*  /  H  (r>Cjo) 

s*r 

Eqs.  (8l)  and  (82)  hold  in  the  interior  of  the  unit  circle.  By  analytic 

continuation  we  can  take  them  to  hold  at  those  points  on  the  unit  circle 

where  the  formulas  exist.  It  is  precisely  at  the  discontinuities  of  the 

loads  that  the  expressions  become  singular. 

On  the  following  pages  we  tabulate  the  closed  form  solutions 

for  the  four  problems  examined.  The  formulas  for  concentrated  load 

agree  with  those  derived  from  the  stress  functions  given  by  Mindlin 
28 

and  Ug  .  The  concentrated  normal  load  disagrees  with  that  deriv'd 
22 

from  f'ichell  ,  but  apparently  a  term  was  left  out  of  hi3  formula. 


DloTHlJUTSD  Uiir  iiORKAL  J'XHiiSJ 


+  + 1 )  C  (r  ,<<  )-4*/rv  oC  4-  CAfi)CCr>0)  ^dCrv/3 

-  3Cr,oC  )(Ccr^oC  -  r  )  -  S(r  >3)( Cxk*3  ~  r )]  j 

=  ^r:([OA+2)rl  +  2(i*:)i,r  +  2Jt] 

+  -£-[Cfc+l)C(r,dCKooa«C  -  r)  -(i.fl)C(r,6)(co^0-r) 

+  5(r,oc)  -  S(r,/J)  0]}  (03) 

=  ITT {  *  -  ^  -  2X  ~  1 

+•  2. [5(r(<<)  +  5(t,3)]  +^r  £K(r,k)-A</n  oc  +■  K(f>B)-A^*v/3j J 

+  2[S(r>oc)  +’'SCr,3)J  —  2 r  Q^(r,^)^*v«c  +  j 


too,  of  -  COO- 
(A+l)r 


-ir-,[Zfr.*c)  ^Ausru  oC  Z. (r,0)^n,jSl  j 


DIJTIUdUTiD  Ui.IT  oiiEAK 


iTHr-oo 


u  =  -dp  (-[^- 2 ^  +  4(*+l)  +  ><A*l r]  Sss0-zs^i 
+  -^-^~^4-l)[C(*’)rf)CoU't<  -  C(  f  »S)  ccra^^J 
+  CrtAr"l)[C(r)oC)-  C  (r,0  )j  —  S(r,<<  )~Aon  aC  +  S(f,0)^»v|^  j 

"  =  4fi>  f T+  ‘+  2 W  iCA  ♦ 

+f  [[(l  +  Or  ->»■'*]  [&(r,<0  +  S(r,£)]  -S(r.rf.)eoorf 
5 C r,  f3 )  Co^Q  +(^+l)[c(r,«<)^»vof  +  C(V,/3)-a<^v<3]]]  j 
°V  =  TTT  ((.A  r  2  -  2  J*  -  +20-  r-'4)[C(rvO  -  C(r<3)] 


+  (^  +  3)  r'1  [K(r,  6)  Cxro/0  -K(Y,<)  cxra<]  (el) 

“(Jk  tl)r"1  |z(r,<0.*t«v<  '2(r,3)^n  0]  +  Ote4l  +2  r'a)[k(ry)  -  k(r,(3)|^ 

+20  +  r‘1)[c(r><)-C(r,^] 

+U+3-2r'l)[K(r,«<)-k(r>6)]-CJb  +  lV1[k(r,-()co^<  -k(r,0)o&0] 
-Cft  +3) r1  [Z(r,t<)^on.oC  -Z(r,0)^on  &]} 


a 


4- t*  i[j?(h*0 Cxro><  +  Z (r,  <9)  Coo-3j 4- (l- If'  ')[Z(^)+ZCr,li'3 
+  r  ^ktr.O -**•*.•<  +k(r,fi)-dim.  fl]  +Zr-I[S(^)+S(r,0')]j 


U2 


CONCENTRATED  UNIT  NORMAL  FORCE 

\1=  4 iry  +  (C-fe-'Ji'it 1  -2Cl  +  0+2(Jbf2)-&g»-]-£|^jy 

+^-[(A4l)C(r,0)c<^e  +  S(r.e-WG]  -(r-  f")  k(r,G)j 

AT  =  ?TT>{-[(3-fc+Z)i-1  -2  +2Ck+2Ekg  r] 

+J.[S(r£)  000,-e-  -Cktl)C(r,Q)  -A4/n.^j  +  (  f  -r-')Z(r,e)] 

°r  =  -(Mr%- 2A  "0(4^7  (( 

4-  2  t  1  [](tr  +  Cxm* 0- )  K(f  ,■ 0-)  -  P(r,-0-)-44/n/&]  J 

<r*=Trr{-l  -(3Ar* 

+  2r-‘  [Z(r,«-)  +  P(r,0)]^n.«-} 

Tr.--M-^rl  +  l)a?fc 

-2f‘[Q(r,«)  0-  +•  Z  (t,-9-)ccra-0-]'| 


u  =  4ir>  lfc*-a>r*  +2a  Ofs^r 

+  I +■  1  )C(r,0-) -d-uw.'9‘  -  S(  cc^O-j 

i-(r-r“1)Zrr,-0-)} 

^  -  4TT>  [  [(3JU2)r*  +  2ftJl+J)  +  2rt+2)^r]  ffj^j 

+  +  l)C(r,d)  co^.0-  -  5(r,©)  —  rj 

+-  (r-i-"1)KCt-1-0-)  -  (r  +  »-  1 )  j 

=Tfr{  (2-fc+l  Vfj77y7 

+  (*♦  r-MZ(r.e-)  t-d-r-1)  P(r.e-) 

■+■  C -fc  +■  3 )  r  k  (r  “(-fcfljf  *  (f.'O’) 

-*-(’-A+l)r  K(r,-0- )^i^, -G-  “(-fe ^3)r  1  2 (r,£) co^'B' j 

T>»=rr  {  (1  fJtr^i^Tr  +  r'1  A  r'1 

+  n-r-J)Q(r,A)  -  r-‘z(r^) 

~A i/n.  A  ■+■  r  k(r  ,-0-)  c/j&Q  ] 


NUMERICAL  REJULTS 


In  a  formal  sense,  the  solution  is  now  completed  and  all 
required  numerical  values  can  be  obtained  from  the  summation  of  the 
series.  In  practice,  however,  this  summation  can  not  be  done  for 
certain  configurations.  In  particular,  if  the  ring  thickness, 

A  =  1  -$  ,  is  small,  convergence  is  so  slow  as  to  be  meaningless. 

We  can  see  this  more  clearly  if  we  look  at  one  of  the  series  expressions 
in  detail. 


S*2 


+  [bt  S  +  (ba  +  4  b  + 

+  [(C,  +-  C*  5  +  C4S  &4SJr”5  1 


(87) 


[>d,  + 


where 


Ds  =e,  +  (t^  +  e3sa)Sxs  +  e4$4*  (88) 

and  where  Qj,  bj  ,  Cj  ,  dj  ,  and  6j  are  constants  dependent  on 
material  properties  and  independent  of  \r  and  5  .  This  series  repre¬ 
sents  the  formula  for  radial  displacement  in  the  annulus  of  a  composite 
disc  of  unit  outer  radius  loaded  by  normal  stress  uniformly  distributed 
over  the  arc  (-♦/z, ♦/*>).  The  formulas  for  other  quantities  are  similar. 


^5 


Terns  of  the  type  C$5  d  f*  ,  for  example,  increase  for  a 

substantial  number  of  terms  before  starting  to  decrease.  For  a  small 
value  of  the  factor  2)  also  increases  with  S  for  many 

terms.  As  a  result  of  these  conditions  convergence  is  very  poor  in 
many  cases.  For  example,  a  case  in  which  Jk  and  4>  each  had  the 
value  .01  was  run  for  five  thousand  terns  with  no  sign  of  convergence. 

A  scheme  was  developed  for  improving  the  convergence  of  the 
series  by  establishing  as  a  first  approximation  the  solution  to  the 
case  for  the  homogeneous  disc.  The  approximation  i3  then  corrected  by 
adding  to  it  the  difference  between  the  corresponding  terms  of  the 
composite  and  the  homogeneous  disc  solutions.  Since  the  terms  of  the 
two  series  approach  one  another  asymptotically  this  offe.s  a  chance 
of  improvement.  Two  forms  of  this  scheme  were  used.  In  one,  the 
material  properties  used  for  the  homogeneous  disc  were  those  of  the 
core  region  of  the  composite  discj  for  the  second  the  properties  used 
were  those  of  the  ring.  In  the  homogeneous  disc,  the  only  elastic 
property  affecting  the  stresses  is  Poisson's  ratio.  This  dependence 
shows  up  in  the  series  in  only  the  term  for  $  =  +l.  Therefore,  regard¬ 
less  of  the  form  used  for  the  first  approximation,  the  speed  of  convergence 
for  the  stress  series  is  the  same.  For  the  displacements,  however, 
the  choice  of  first  approximation  does  make  a  difference.  It  was  found 
that  about  half  as  many  terms  are  required  if  the  ring  properties  are  used. 
The  bulk  of  the  computations,  however,  had  been  performed  using  the 
properties  of  the  core.  These  computations  were  spot  checked  using  the 


more  rapidly  convergent  technique 
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Numerical  values  have  been  outained  for  selected  points 
on  tne  outer  surface  and  along  tne  interface  for  the  case  of  unit 
normal  load  acting  on  a  composite  disc  in  olane  stress.  Two  config¬ 
urations  of  material  constants  nave  been  used,  corresponding  respec¬ 
tively  to  a  layer  of  gold  or  alumina  on  a  base  of  steel.  The  proper¬ 
ties  of  tnese  materials  are  shown  below  in  Table  I. 


Table  I.  Material  Properties 


Materi al 

A*(psi) 

h) 

A 

y/y 

Otcel 

1.153(10) ' 

0.300 

1.86 

1.000 

r.ola 

(i 

U. 014(10) 

0.1*20 

2, 1*5 

0.3^81 

Alumina 

p.086(10)7 

0.2U6 

1.31' 

2.1*18 

Values  are  calculated  for  a  dimensionless  configuration  of 
unit  radius,  unit  shear  modulus  in  the  base  region,  and  unit  normal 
stress.  *'e  can  convert  results  to  a  system  with  dimensionalized 
quantities.  Let  and  d-(r.e-)  be  the  displacements 

a 

and  stresses  in  a  disc  of  radius  p  ,  witn  a  snear  modulus  ^u  in  the 
core,  rfe  can  convert  the  dimensionless  results  to  values  for  a 
corresponding  dimensionalized  disc  by  formulas  such  as 

u  (r.O)  =  -p-  u(  r/ p 

(  :?) 

0(j  ,9-)  =  0-  (r/e,0) 

Table  II  lists  some  results  for  the  disc  whose  layer 
thickness  J\.  and  load  width  4*  both  are  equal  to  0.01.  Values  in  the 
vicinity  of  tne  .oad,  both  on  the  surface  anu  nt  tne  interface,  are 


given.  A  jump  discontinuity  in  the  hoon  stress  at  the  interface  is 

A  I 

indicated  by  listing  as  the  value  in  the  ring  and  as  CF^. 

the  value  in  tne  core.  A  list  is  shown  for  botn  a  gold  (soft)  and 
alumina  (hard)  layer  on  a  steel  core.  Figure  2  shows  graphically 


these  and  some  other  results 


« I 
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Figure  2 
Part  D 


Behavior  of  a  Composite  Disc 
T..Aat  the  Interface. 
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o.oo8 


0.010 


Figure  2 
Part  C 


Behavior  of  a  Composite  Disc 
O ’  at  the  Interface. 


Figure  2.  Behavior  of  a  Composite  Disc 
Part  E.  at  the  Interface  in  the  Base. 
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Figure  2. 
Part  F, 


Behavior  of  a  Composite  Disc 
0[^at  the  Interface  in  the  Ring. 

Tr 


Table  II 


DISPLACEMEXTS  AND  STRESSES* 
Part  A.  Gold  Layer  on  Steel 


100  Surface  Tnterface 


X 

* 

100  u 

100  u 

gy 

TV* 

_  A 

5±L 

0.00 

1.833 

0.797 

1.096 

0.612 

0.000 

0.219 

0.083 

0.05 

1.831 

0.798 

1.096 

0.610 

-0.023 

0.218 

0.084 

0.10 

1.823 

0.799 

1.094 

0.603 

-0.045 

0.217 

0.086 

0.15 

1.809 

0.802 

1.091 

0.594 

-0.066 

0.215 

0.089 

0.20 

1.790 

0.805 

1.087 

0.580 

-0.087 

0.212 

0.092 

0.25 

1.763 

0.810 

1.081 

0.563 

-0.106 

0.209 

0.097 

0.30 

1.730 

0.815 

1.075 

0.543 

-0.123 

0.205 

0.103 

0.35 

1.686 

0.821 

1.068 

0.519 

-0.139 

0.200 

0.109 

0.40 

1.631 

0.828 

1.060 

0.494 

-0.153 

0.195 

0.116 

0.45 

1.560 

0.835 

1.051 

0.467 

-0.164 

0.189 

0.123 

0.50 

1.441 

m 

1.042 

0.438 

-0.173 

0.183 

0.130 

0.55 

1.321 

-0.148 

1.032 

0.408 

-0.180 

0.177 

0.137 

0.60 

1.247 

-0.140 

1.021 

0.378 

-0.184 

0.170 

0.143 

0.65 

1.188 

-0.131 

1.011 

0.349 

-0.187 

0.163 

0.149 

0.70 

1.138 

-0.122 

1.000 

0.319 

-0.187 

0.156 

0.155 

0.75 

1.097 

-0.113 

0.988 

0.291 

-0.185 

0.150 

0.159 

0.80 

1.060 

-0.103 

0.977 

0.264 

-0.182 

0.143 

0.163 

0.85 

1.029 

-0.094 

0.965 

0.238 

-0.177 

0.136 

0.166 

0.90 

1.001 

-0.086 

0.954 

0.214 

-0.172 

0.129 

0.168 

0.95 

0.976 

-0.077 

0.943 

0.191 

-0.165 

0.123 

0.169 

1.00 

0.954 

-0.069 

0.932 

0.170 

-0.158 

0.116 

0.169 

1.25 

0.871 

-0.033 

0,860 

0.093 

-0.118 

0.089 

0.158 

1.50 

0.817 

-0.009 

0.834 

0.050 

-0.084 

0.068 

0.138 

1.75 

0.778 

0.004 

0.795 

0.028 

-0.058 

0.053 

0.116 

2.00 

0.747 

0.011 

0.762 

0.017 

-0.041 

0.042 

0.097 

2.25 

0.721 

0.013 

0.732 

0.011 

-0.030 

0.034 

0.081 

2.50 

0.697 

0.012 

0.706 

0.008 

-0.022 

0.028 

0.068 

2.75 

0.676 

0.010 

0.683 

0.006 

-0.017 

0.024 

0.058 

3.00 

0.656 

ts.008 

0.661 

0.005 

-0.013 

0.020 

0.050 

3.25 

0.637 

0.006 

0.642 

0.004 

-0.011 

0.018 

0.043 

3.50 

0.620 

0.005 

0.623 

0.003 

-0.009 

0.015 

0.038 

3.75 

0.603 

0.003 

0.606 

0.003 

-0.007 

0.013 

0.033 

4.00 

0.586 

0.002 

0.591 

0.002 

-0.006 

0.012 

0.029 

4.25 

0.574 

0.001 

0.576 

0.002 

-0.005 

0.010 

0.026 

4.50 

0.560 

0.001 

0.562 

0.001 

-0.004 

0.009 

0.023 

4.75 

0.547 

0.000 

0.549 

0.001 

-0.003 

0.008 

0.021 

5.00 

0.535 

0.000 

0.536 

0.001 

-0.003 

0.007 

0.018 

*in  a  composite  disc  of  unit  radius  loaded  by  a  unit  tensile  radial  stress 
over  the  arc  (-0.005,  0.005)  and  by  a  concentrated  reaction  at  the  center. 
Layer  thickness  is  0.01. 
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Table  II 

DISPLACEME1TS  AID  STRESSES* 
Pert  B.  Alumina  Layer  on  Steel 


100  Surface  Interface 


X 

4 

100  u 

100  U 

°~r 

tr. 

__  A 

I 

Ojr 

0.00 

1.151 

1.271 

1.020 

0.485 

0.000 

-0.227 

0.006 

0.05 

1.150 

1.271 

1.020 

0.483 

-0.017 

-0.225 

0.006 

0.10 

1.148 

1.269 

1.018 

0.480 

-0.034 

-0.220 

0.008 

0.15 

1.145 

1.266 

1.017 

0.474 

-0.051 

-0.211 

0.011 

0.20 

1.139 

1.262 

1.014 

0.465 

-0.068 

-0.199 

0.015 

0.25 

1.133 

1.257 

1.010 

0.455 

-0.082 

-0.184 

0.019 

0.30 

1.124 

1.250 

1.006 

0.442 

-0.096 

— 0.166 

0.025 

0.35 

1.114 

1.243 

1,002 

0.428 

-0.109 

-0.146 

0.031 

0.40 

1.101 

1.236 

0.996 

0.412 

-0.120 

-0,124 

0.038 

0.U5 

1.084 

1.227 

0.990 

0.395 

-0.130 

-0.100 

0.045 

0.50 

1.059 

• 

0.984 

0.377 

-0.138 

-0.075 

0.052 

0.55 

1.033 

0.208 

0.977 

0.358 

-0.145 

-0.050 

0.060 

o.6o 

1.015 

0.198 

0.970 

0.339 

-0.150 

-0.024 

0.067 

0.65 

0.999 

0.187 

0.963 

0.320 

-0.154 

0.001 

0.075 

0.T0 

0.985 

0.176 

0.955 

0.301 

-0.157 

0.025 

0.08? 

0.75 

0.972 

0.165 

0.947 

0.282 

-0.158 

0.048 

0.088 

0.60 

0.960 

0.154 

0.939 

0.263 

-0.158 

0.070 

0.094 

0.05 

0.949 

0.142 

0.931 

0.246 

-0.157 

0.090 

0.099 

0.90 

0.938 

0.131 

0.922 

0.229 

-0.155 

0.108 

0.104 

0.95 

0.927 

0.121 

0.914 

0.213 

-0.152 

0.124 

0.108 

1.00 

0.917 

0.110 

0.906 

0.198 

-0.149 

0.139 

0.1U 

1.25 

0.871 

0.063 

0.865 

0.135 

-0.127 

0.186 

0.119 

1.50 

0.830 

0.028 

0.827 

0.093 

-0.104 

0.201 

0.116 

1.75 

0.794 

0.005 

0.792 

0.064 

-0.084 

0.197 

0.107 

2.00 

0.761 

-0.010 

0.761 

0.044 

-0.068 

0.185 

0.097 

2.25 

0.731 

-0.017 

0.732 

0.031 

-0.055 

0.169 

0.087 

2.50 

0.705 

-0.020 

0.706 

0.021 

-0.044 

0.153 

0.077 

2.75 

0.681 

-0.020 

0.682 

0.015 

-0.036 

0.137 

0.068 

3.00 

0.659 

-0.019 

0.660 

0.010 

-0.030 

0.122 

0.060 

3.25 

0.639 

-0.017 

0.640 

0.007 

-0.024 

0.109 

0.053 

3.50 

0.621 

-0.014 

0.622 

0.005 

-0.020 

0.097 

0.047 

3.75 

0.604 

-0.012 

0.605 

0.003 

-0.017 

0.087 

0.042 

4.00 

0.588 

-0.010 

0.589 

0.002 

-0.015 

0.078 

0.037 

4.25 

0.573 

-0.008 

0.574 

0.002 

-0.012 

0.070 

0.033 

4.50 

0.559 

-0.006 

0.560 

0.001 

-0.011 

0.063 

0.030 

4.75 

0.546 

-0.005 

0.547 

0.001 

-0.009 

0.058 

0.027 

5.00 

0.534 

-0.004 

0.534 

0.001 

— 0.008 

0.052 

0.025 

*in  a  composite  disc  of  unit  radius  loaded  by  a  unit  tensile  radial  stress 
over  the  arc  (-0,005,0.005)  and  by  a  concentrated  reaction  at  the  center. 
Layer  thickness  is  0.01. 
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PART  V. 

DEFORMATION  COEFFICIENT  OF  FRICTION 

If  a  load  is  applied  to  a  portion  of  the  surface  of  a 

body  there  vill  be  an  associated  deformation.  There  will  then  be 
stored  in  the  body  strain  energy  which  is  equal  to  the  work  done  by 
the  load  in  deforming  the  body.  We  shall  make  certain  assumptions 
regarding  the  behavior  of  this  strain  energy  when  a  load  is  trans¬ 
lated  an  infinitesimal  distance  on  the  surface.  From  these  assump¬ 
tions  we  vill  determine  a  friction  force  which  leads  to  the  idea  of 
a  deformation  coefficient  of  friction.  It  must  be  emphasized  that 
this  deformation  concept  is  not  suggested  as  an  alternative  to  the 
more  commonly  suggested  mechanisms  such  ns  rupture  of  asperities  and 
sticking.  Rather  it  complements  them  as  an  additional  possible 
mechanism. 

Consider  the  case  of  a  constant  normal  stress  acting  on  a 
surface  element  of  the  disc.  A  neighboring  state  of  loading  is 
defined  as  a  condition  in  which  the  load  has  been  translated  circum¬ 
ferentially  an  incremental  distance,  say  Sx  .  Clearly  the  amount 
of  strain  energy  in  the  disc  is  unchanged  as  the  system  goes  from 
one  state  of  loading  to  a  neighboring  state.  In  a  local  sense, 
however,  there  has  been  a  change.  The  radial  displacement  of  each 
point  on  the  surface  has  changed  to  that  of  its  neighbor.  Part  of 
the  disc  has  increased  in  strain  energy  while  some  has  decreased. 
Specifically,  a  certain  amount  of  work  is  done  by  the  normal  stress 
acting  on  that  part  of  the  loaded  segment  in  which  the  magnitude  of 


the  deformation  increases.  An  equal  amount  of  work  is  done  by  the 
deformed  body  over  that  part  in  which  the  deformation  decreases. 

The  basis  for  the  postulate  of  a  deformation  coefficient  of 
friction  stems  from  the  assumption  that  the  strain  energy  added  to 
the  system  is  supplied  by  the  action  of  a  circumferential  force, 
while  the  strain  energy  released  by  the  system  is  dissipated  in  the 
form  of  heat,  noise,  or  other  non-re vers able  action. 

Although  the  model  does  not  distinguish  between  sliding  and 
rolling,  the  phenomenon  is  somewhat  easier  to  visualize  if  we  think 
of  a  rigid  roller  indenting  a  surface  as  it  rolls.  Work  is  required 
to  deform  the  surface  ahead  of  the  roller.  The  energy  released  by 
the  restoration  of  the  initial  configuration  behind  the  roller  is 
assumed  lost.  Simkins  suggests  that  interaction  of  molecules  of  the 
mating  surfaces  during  sliding  causes  an  additional  dissipative 
effect  not  accounted  for  in  this  model.  This  interaction  would  be 
missing  during  rolling,  however.  He  suggests,  therefore,  that  the 
model  is  more  realistic  when  thought  of  in  terms  of  rolling  contact. 

To  calculate  the  deformation  coefficient  we  shall  consider 
the  elastic  deformation  of  the  surface  due  to  a  normal  stress  0  • 

With  reference  to  Figure  3*  divide  the  loaded  half-interval  into  rt 
segments  of  width  &X  ,  where  &X  is  the  amount  of  the  circum¬ 
ferential  displacement.  Associated  with  this  displacement,  the 
radial  deformation  of  point  X  i  becomes  U  .  The  work  required 
to  deform  the  loaded  surface  in  the  l  th  interval  is 


"Thomas  E.  Simkins:  Private  communication 
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The  expression  for  the  last  interval  is  one  half  of  this.  Summing 
over  the  pertinent  half-range  and  simplifying  gives 

SW  Awl  =  9  i  x  [•“« Y u»  -  ](!>1) 

Since  lx  is  small  and  1\  is  large  ,  ve  may  write 

S  W  =9  lx  (u0  -  '9*> 

The  force  P  required  to  effect  the  translation  must,  by 
hypothesis,  provide  the  energy  $W.  Therefore 

P  =  SW/ix  =  9  (U0  -  un)  (93) 

In  a  disc  of  radius  p  loaded  by  a  uniform  normal  stress 
0  over  an  arc  of  width  4>  ,  the  total  force  F  exerted  by  the  load 
is 

F=9 

The  coefficient  of  friction  F  is  defined  similar  to  the  usual  form 

where  the  absolute  value  must  be  taken  since  the  model  does  not 
distinguish  between  tensile  and  compressive  forces.  Since  the  defor¬ 
mation  is  directly  proportional  to  the  intensity  0  of  the  stress, 
the  frictjcn  coefficient  is  a  linear  function  of  the  normal  force 


applied  over  a  constant  element  of  arc. 

To  illustrate  the  computation ,  let  us  again  use  the  caret 
to  distinguish  dimensional  quantities.  Then 


u  =  (g  p/yu  )  u(l,G) 

and 


(96) 


(97) 


Applying  the  data  for  gold  on  steel  as  shown  in  Table  II ,  we  obtain 

T  =  3-  +  <10)'*|g|  (98) 

The  friction  coefficient  due  to  the  deformation  is  therefore 
quite  small.  For  loads  of  moderate  intensity,  the  effect  is  negligible 
compared  to  other  friction  mechanisms. 

When  the  loaded  arc  width  is  small  the  friction  coefficient 
is  essentially  independent  of  for  a  fixed  ratio  Jk/<P  .  Results  will 
be  given  in  texms  of  the  DFC,  that  is  the  coefficient  which  multiplies 
the  stress  to  obtain  the  conventional  coefficient  of  friction.  Shown  in 
Figure  4  is  the  behavior  of  the  DFC  for  ^  =O.Ol#  for  both  cases,  gold 
on  steel  and  alumina  on  steel.  Table  III  lists  some  values  of  the 
normalized  friction  coefficient  for  these  two  cues. 
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Table  III 

DEFORMATION  COEFFICIENT  OF  FRICTION 
Part  A.  Gold  Layer  on  Steel  Base 


Q  _ 

DFC  X  (10)  (in  /lb)  for  prescribed^ 


A/<P 


0.0  (All  Base) 

1.23 

0.005 

1.24 

0.01 

1.25 

0.02 

1.27 

0.05 

1.33 

0.1 

1.1*3 

0.2 

1.65 

0.5 

2.64 

1.0 

3.31 

2.0 

- 

5.0 

- 

10. 

- 

20. 

- 

50. 

- 

100. 

- 

200. 

- 

VJ1 

o 

o 

• 

- 

(All  Layer) 

3.31 

0 


0.1 

0.01 

0.001 

1.45 

1.47 

1.47 

1.46 

1.49 

- 

1.47 

1.49 

- 

1.49 

1.51 

- 

1.56 

1.58 

- 

1.67 

1.69 

1.70 

1.90 

1.92 

1.93 

2.70 

2.71 

2.71 

3.40 

3.40 

3.40 

3.71 

3.73 

3.73 

3.81 

3.84 

3.85 

3.82 

3.86 

3.86 

- 

3.86 

3.87 

- 

3.86 

3.87 

- 

3.86 

3.87 

- 

- 

3.87 

- 

- 

3.87 

3.82 

3.86 

3.87 
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Table  III 

DEFORMATION  COEFFICIENT  OF  FRICTION 
Part  B.  Alumina  Layer  on  Steel  Base 

DFC  XdO)8  (in2/lb)  for  prescribed 


<P 


1/<P 

1.0 

0.1 

0.01 

0.001 

0.0  (All  Base) 

1.23 

1.1*5 

1.1*7 

1.1*7 

0.005 

1.22 

l.l*J* 

- 

- 

0.01 

1.22 

1.1*1* 

1.1*6 

- 

0.02 

1.21 

1.1*2 

1.1*5 

- 

0.05 

1.18 

1.1*0 

1.1*2 

- 

0.1 

1.15 

1.36 

1.38 

- 

0.2 

1.06 

1.26 

1.28 

m 

0.5 

0.72 

0.97 

0.99 

0.99 

1.0 

0.56 

0.78 

0.80 

0.80 

2.0 

- 

0.70 

0.72 

0.72 

5.0 

m 

0.68 

0.69 

0.70 

10. 

- 

0.68 

0.69 

0.69 

20. 

- 

- 

0.69 

0.69 

50. 

- 

- 

0.69 

0.69 

100. 

- 

- 

0.69 

0.69 

200. 

- 

- 

- 

0.69 

500. 

- 

- 

- 

0.69 

(All  Layer) 

0.56 

0.68 

0.69 

0.69 

ft'  *?' 
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PART  VI. 

CONCLUSIONS 

The  finite  Fourier  transform  has  been  used  to  develop  a 
field  solution  for  a  composite  disc  in  a  at  ate  of  plane  elastic 
equilibrium.  The  method  has  the  advantage  of  being  explicit  compared 
to  more  familiar  methods.  No  prior  knowledge  of  the  solution  is 
required.  On  the  other  hand  it  has  the  disadvantage  of  being  somewhat 
cumbersome,  and  it  does  not  make  use  of  the  accumulated  collectiona  of 
known  solutiona  to  the  biharmonic  problem. 

The  properties  of  the  finite  Fourier  transform  cause  the 
method  to  have  merit  for  circular  regiona  with  all  boundaiys  concentric 
circles.  The  problem  examined  in  this  investigation  was  a  fundamental 
elasticity  problem  of  the  flrat  kind.  If  displacements  were  prescribed 
Instead  of  streases  the  technique  should  be  equally  useful.  In  the 
event  of  mixed  boundary  conditions,  however,  the  method  may  lead  to 
dlfflcultlea  asaociated  with  transforming  the  boundary  conditions. 

In  the  preaent  problem  the  transform  of  the  pair  of  partial 
differential  equations  in  two  independent  displacement  variables 
reduces  to  a  pair  of  ordinary  differential  equations.  These  are 
solved  explicitly.  The  inversion  of  the  transformed  solution  is 
obtained  by  aumaing  the  series  of  solutions  to  the  transformed  problem. 

The  numerical  evaluation  of  the  sum  was  accomplished  with 
some  difficulty.  In  most  of  the  cases  of  interest,  direct  evaluation 
could  not  be  performed.  Convergence  was  obtained,  however,  by  using  the 
closed  fora  solution  for  a  homogeneous  disc  as  a  first  approximation. 
Even  then,  convergence  is  slow.  Because  an  increasingly  large  musber 
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of  terms  is  required  as  the  loaded  sector  becomes  small,  there  may  be 
an  accumulated  rounding  error  for  some  results.  Double  precision 
computation  could  eliminate  some  of  these.  If  results  are  required 
only  in  the  vicinity  of  the  load,  the  solution  to  the  problem  of  a 
semi-infinite  plate  vould  be  adequate  for  the  small  loaded  sector 
vhere  curvature  effects  are  negligible.  It  was  found  that  for 
*$*  $  0.  1  one  could  essentially  disregard  curvature  effects. 

A  postulated  mechanism  for  friction  on  a  smooth  surface  was 
developed.  Friction  results  from  assumptions  about  the  dissipation 
of  strain  energy  when  a  deformed  body  is  transformed  to  a  neighboring 
deformation  state.  The  friction  coefficient  found  by  this  method  is 
very  small}  its  presence  vould  be  obscured  by  more  significant  factors. 
However,  this  frictional  contribution  could  be  a  limitation  on  the 
complete  elimination  of  friction  between  solid  surfaces.  To  the 
extent  that  this  theory  ii  valid,  bodies  with  a  high  elastic  modulus 
vill  exhibit  less  friction  than  one  with  a  low  modulus  due  to  the 
smaller  amount  of  strain  energy  associated  vith  a  given  load. 
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